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Abstract 


The following applications of drops and bubbles in materials science are reviewed s 
Formation of extended p-n junctions in semi-conductors by drop migration; mechanisms and 
morphologies of migrating drops and bubbles in solids; nucleatlon and corrections to the 
Volmer-Weber equations; bubble shrinkage in the processing of glass; the formation of glass 
microshells as laser-fusion targets; and radiation-induced voids in nuclear reactors. These 
examples show that knowledge of the behavior of drops and bubbles is .' iportant in diverse 
areas of materials processing and properties, and deserves further fundamental study. 

Introduction 


Semiconducting junctions can be made in complicated shapes by migration of liquid bubbles 
in a solid semiconductor. Nucleatlon of solid phases must be controlled in processing solids 
and regulating their properties; the theory of nucleatlon was first derived for the forma- 
tion of droplets in a vapor. Melting of clear and strong glass requires that bubbles be 
removed from viscous molten glass, and tiny microshells of glass as laser fusion targets are 
formed by blov/ing a bubble inside a drop of molten glass. These examples show the impor- 
tance of drops and bubbles in processing solid materials. Radiation can form bubbles in 
reactor materials, weakening them. In each of these areas there have been recent advances 
or controversies, which will be described in this paper. These examples are not intended 
to provide a comprehensive review of drops and bubbles in Materials Science, but merely a 
sample of applications and opoortunities. 

Semiconductor processing and migration of drops 

Commercial processes for making p-n junctions in semiconductors by migration of li?^uid 
drops, wires, or sheets have been developed by T. R. Anthony and H. E. Cline of Ge'' 

Electric. Molten metallic drops or wires move through a solid semiconductor ' lermal 

or electrical gradient, leaving resolidified doped material in their wake. Dot rlpes 

of a low-melting metal that is a suitable dopant, such as aluminum, are deposit the 

surface of the semiconductor, which is then heated and subjected to a temperatitre . adient 
or electrical field. In practice there are a number of problems in controlling the shape 
of the liquid region and its path that require careful control of the temperature and the 
gradients. Thus knowledge of the morphology and kinetics of the liquid drops and wires is 
important for successful production of a satisfactory device. Thermal gradients proved 
be difficult to control, so electromigration was developed as an alternative.^ 

A lattice of columnar p-n junctions through a semiconductor wafer can be made by thermal 
migration of aliiminum drops through silicon.' This "deep-diode array" can be used as an 
imaging target for x-rays and infrared radiation. 

A liquid drop in a gradient of driving force in a solid moves by dissolution and repreci- 
pltatlon of solid in the liquid. Driving gradients can be in temperature, mechanical force 
(acceleratlonal field) or electrical potential. Anthony and Cline found that high force 
gradients were required to migrate brine drops in solid KCi,^ suggesting that an imprac- 
tlcally high sedimentation field would be required to move liquid drops in semiconductors. 

The rate of droplet motion can be controlled either by the rate of diffusion oi solute 
in the liquid or by the rate of dissolution of solid (interface control) . This dissolution 
rate is usually proportional to the solubility of the solid in the liquid drop, and can be 
different for different crystallographic planes in the solid. At high temperatures con- 
vective flow in the drop can lead to rates higher than expected from the static diffusion 
coefficient of solute in the liquid.^ For large drops diffusion in the liquid is more 
likely to control the rate of droplet migration, whereas for small drops interface rates 
become important. A similar influence of crystal size is found for the rate of crystalli- 
zation from solution. For large crystals diffusion in the liquid controls the rate of 
growth, whereas for small crystals interface control is usual. ^ 

At equilibrium the shape of a liquid drop in a solid is cok.crolled by the int&rfacial 
energies of different crystallographic planes of the solid in contact with the liquid. The 
shape that minimizes the total surface free energy is the equilibrium one. In principle 
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this ahape can be calculated from knowledge of the surface energy as a function of crystal- 
lographic orientation of a surface facet by the "Wulff" construction; In practice this 
function is hardly ever known in much detail, and the relative values of Interfaclal 
energies are often deduced from equilibrium shapes of crystals. One would expect that low 
index and closely-packed crystallographic planes should have the lowest energies, so that 
these planes should, and usually do, predominate in equilibrium shapes of crystals. 

Brine droplets in KCt had mainly (100) planes as faces after seven years at room tempera 
ture, with edge roundings and some contribution from (110) faces. ^ Silicon ex^ses mainly 
(111) planes in contact with liquid, with some contribution from (100) planes.' ° When a 
droplet migrates in a force field, its shape can be distorted by nonuniformities in the 
field and in the rates of interfacial dissolution.'” 

The rate of migration of gac-filled bubbles in KCe in a temperature gradient was deter- 
mined by the rate of evaporation of solid and by diffusion of KCt molecules in the bubble, 
in much the way as the liquid drops." The bubbles were bounded by (100) planes, but 
stretched out in the direction of migrat >n instead of remaining cubic." The shape was 
determined by the interfacial rates of evaporation and deposition rather than by the equili 
brium condition from surface energies. A trail of fine bubbles followed in the wake of the 
migrating bubble. 


Nucleatlon 


The transformation of one phase into another almost always requires formation of a mole- 
cular cluster of the new phase in the old, after which the new phase giows at the expense 
of the old. This nucleation and growth mechanism is valid for pure materials and usually 
for mixtures; in some mixtures a uniform transformation called soinodal decomposition is 
possible, but it will not be considered here. When a cluster or nucleus foms the resul- 
tant surface energy can exceed the volume free energy difference driving transformation, 
so that a nucleation barrier exists even when a large sample of the new phase is stable 
with respect to the old. This situation is illustrated in Figure 1, where the chemical 



Figure 1. The chemical potential of ethanol drop at -9*C, as a function of drop radius. 

potential in Vntolecules of a drop of ethanol is compared to the chemical pote.itial of 
supersaturated ethanol vapor uy and the chemical potential of liquid ethanol ue, all at 
-9*c. When the chemical potential of a drop equal^i that of the surrounding vapor it is 
said to be of critical size (radius r*) . Smaller drops are unstable and decompose; larger 
drops grow. 

Volmer and Weber derived the number of critical nuclei as a function of superjaturation 
and temperature from the fluctuation theory of Einstein" and the calculation of Gibbs" 
of the reversible work to form a sphere of the critical size. The result is 

N* - N exp(-W/kt) (1 
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where N* is the number of critical nuclei per unit volume, N is the number of 
cules per unit volume, W is the reversible %iork to form the critical nucleus, 
Boltzmann's constant and T the temperature. Then the rate of nucleation I is 
of the number of collisions of vapor molecules with the surface of a critical 
N*: 


vapor mole- 
k is 

the product 
nucleus wxth 


I = 2 A* N* (2) 

where Z is the collision frequency of molecules with a surface and A* is the surface area of 
a critical nucleus. For a spherical nucleus: 

W = 16nYV3(aP)^ (3) 

where y is the surface energy of the liquid and aP is the difference in pressure bettreen the 
inside and outside of the critical nucleus. 

The equations of Volirer and Weber agreed well with experimental data in rain formation of 
various liquids in a cloud chanfcer^^ and with more recent results in the diffusion cloud 
chamber, especially in the calculated and measured values of the surface energies. Never- 
theless a number of authors have suggested that the equations of Volmer and Weber need to 
be modified by correction factors that vary from a factor of 10^^ in the nucleation rate to 
a negligible correction. These treatments involve questionable assumptions of equili- 
brium in chemical deposition reactions to form the critical nucleus or between embryos not 
of the critical size and the surrounding vapor. The equations of /olmer and Weber agree 
well with experiment, and there appears to be no need to correct them. 

In a recent paper Carlon has found that even undersaturated water vapor and moist air 
contain multimolecular clusters of water. Thus the initial condition of a vapor in a 
nucleation experiment may not be just individual molecules, but could contain an unknown 
distribution of clusters. Such a distribution would further con^licate nucleation theories 
based on a series of condensation reactions, but should not influence the equilibrium 
fluctuation method of Volmer and Weber. Carlon also found that these clusters of water 
molecules lead to infrared absorption at wave lengths different from those 'or molecular 
water. 


Bubbles in the processing of glass 

Pining, or the removal of bubbles from a glass melt, is one of the major technological 
problems in glass melting. It is usually solved by holding the glass for some time at a 
temperature somewhat below the hiahest melting temperature, and by adding certain minor 
constituents to the original glass batch. The mechanisms by which these additions a^d 
bubble removal are still somewhat uncertain, although much progress in understanding has 
been made in recent years. 

Bubbles can be removed from a melt by either of two ways. They cai» rise to the surface 
or the gas in them can dissolve in the glass. The rate of rise is given by the following 
equation: 

^ (A) 

dt 9n 

where c is the density of the glass, g is the gravitational constant, R is the bubble radius, 
and n is the viscosity of the glass. For a viscosity of 100 P, typical for melting tempera- 
tures, the rate of rise of bubbles 0.1 mm in diameter is about 10 cm/day, which is too small 
to eliminate them from a normal glass furnace. Thus small bubbles can be removed from glass 
melts only by dissolution of their gas into the glass melt, although larger bubbles can rise 
to the surface. 

Arsenic oxide is a common fining agent added to glass to help remove bubbles. For many 
years it was thought that the arsenic released oxygen at glass melting temperatures, which 
"swept out” the bubbles in the glass. However, the calculation above shows that such a 
mechanism would not eliminate small bubbles, and the elegant experiments of Greene and 
co-workers showed that arsenic enhances dissolution of oxygen bubbles in glass. Thus 
the importance of arsenic and antimony oxide additions to the glass is to aid in removal of 
fine bubbles, rather than to generate more gas. 

Studies of bubble shrinkage in glass can help in understanding the kinetics of fining, 
the rates of diffusion of gases in molten glass, and the theoretical treatment of bubble 
dissolution. Greene and his co-workers observed the shrinkage of oxygen and sulfur dioxide 
bubbles in a munber of silicate glasses between 1000*C and 1300’C.^^~^^ They sealed a small 
amount of gas into a hole between two glass discs, which were then heated in a rotating 
silica cell. During the first few minutes of heating the temperature of the bubble changed. 
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so that reliable measurements were not possible. This technique also often introduced a 
small amount of foreiqn qas into the bubble. To avoid these difficulties Greene and Davis 
developed a technique in which the bubble was observed just after it was blown into the 
molten glass, and studied oxygen, nitrogen, and water bubbles in molten boric oxide.-** 
Brown and Doremus-'*' also used this technique to study the shrinkage of air oxygen, and 
nitrogen bubbles in molten boric oxide; an example of their work is shown in Figure 2. 



Figure 2 . Diameter of an air bubble in molten B20^ at b50®C as a function of time. Circles 
data; line from eg. 5 with D = 2.3(10)~‘ cm^/sec and B » 0.058. 

The mathematical treatment of the growth or shrinkage of a bubble requires solution of a 
non-linear oartial differential equation, and a general analytical solution has not been 
found. However, several limiting cases can be solved satisfactorily. If the gas dissolving 
in the glass has a low concentration (dilute solution in the glass) , a relatively simple 
relation for the radius of the bubble R as a function of time t is f ound : ^ 

2 2 

R = 2DBt(l + -— ^) (5) 

o 


where Rq is the initial bubble radius, D is the diffusion coefficient of the dissolving gas 
in the glass, and B is a concentration factor B = (c^-c„.) /Cg, wliere the subscripts denote 
concentrations of gas: s in the bubble, i dissolved in the glass at the bubble-glass inter- 
face, and dissolved in the glass far from the bubble. Eq. 5 fits data for the dissolution 
of air i*. molten boric oxide as shown in Figure 2. 

If the amount of qas dissolved in the glass is not small, (B 1), the solution becomes 
more complicated, because the flow of molten glass as the bubble shrinks leads to a "convec- 
tion" term in the differential equation.-' True convection in the glass melt is not sig- 
nificant because of the high viscosity of the glass. A perturbation solution of this 
expanded differential equation for bubble growth or shrinkage has been given by Duda and 
Vrentas;- '* the first term in their solution is eq. 5. 

Further complications arise if more than one kind of qas is dissolving in the 
glass, if the glass sample is not large compared to the bubble,-'^ and if there is 
a surface chemical reaction of the gas dissolving in the glass.’* 

Hollow glass shells filled with deuterium and tritium can serve as laser-fusion targets. 
These shells are tynically about 100 wm in diameter with a wall thickness of a few microns, 
and can be made by dropping oowder of a metal-organic gel through a furnace tower at temper- 
atures of from 1200*C to ISOO’C.*- The organic material and residual water serve to blow 
out the glass into a spherical shell. The size and quality of the shells are dependent on 
the glass composition, the organic content of the gels, and the furnace conditions. 
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Swelling of materials in fast reactors 

The swelling of radiation-induced voids in metals and alloys in nuclear reactors can 
lead to serious degradation of these materials, and may become a problem in fast breeder 
reactors. The rate of stielling is strongly dependent upon the composition of the irradiated 
alloy; minor additions ("swelling inhibitors") can reduce the amount of swelling by up to 
two orders of magnitude. Nucleation of the voids may be a critical step in their forma- 
tion.^** Bubbles of gas in irradiated metals can grow by migration and coalescence. The 
rate of migration of small bubbles is controlled by surface diffusion, and in larger 
bubbles by ledge nucleation, in niobium alloys and zirconium. 
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